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Many optimization problems can be naturally represented as (hyper) graphs, where vertices corre¬ 
spond to variables and edges to tasks, whose cost depends on the values of the adjacent variables. 
Capitalizing on the structure of the graph, suitable dynamic programming strategies can select certain 
orders of evaluation of the variables which guarantee to reach both an optimal solution and a minimal 
size of the tables computed in the optimization process. In this paper we introduce a simple algebraic 
specification with parallel composition and restriction whose terms up to structural axioms are the 
graphs mentioned above. In addition, free (unrestricted) vertices are labelled with variables, and the 
specification includes operations of name permutation with finite support. We show a correspon¬ 
dence between the well-known tree decompositions of graphs and our terms. If an axiom of scope 
extension is dropped, several {hierarchical) terms actually correspond to the same graph. A suitable 
graphical structure can be found, corresponding to every hierarchical term. Evaluating such a graph¬ 
ical structure in some target algebra yields a dynamic programming strategy. If the target algebra 
satisfies the scope extension axiom, then the result does not depend on the particular structure, but 
only on the original graph. We apply our approach to the parking optimization problem developed 
in the ASCENS e-mobility case study, in collaboration with Volkswagen. Dynamic programming 
evaluations are particularly interesting for autonomic systems, where actual behavior often consists 
of propagating local knowledge to obtain global knowledge and getting it back for local decisions. 


1 Introduction 

Many optimization problems can naturally be represented as hypergraphs, where each hyperedge is 
an atomic subproblem and is attached to vertices corresponding to the problem’s variables. However, 
these hypergraphs often lack an algebraic structure, which would allow for the recursive resolution of 
problems, and are not able to represent the secondary optimization problem, that is: finding an optimal 
variable elimination strategy. The order in which variables are eliminated may dramatically affect the 
computation cost of solutions. 

In this paper we introduce a simple algebraic specification for representing optimization problems. 
It is similar to a process calculus, based on nominal structures (namely permutation algebras, see e.g. 
milTdl). Optimization problems are represented as terms over the variables of the problem, consisting 
of the parallel composition of subproblems. A key feature is the representation of variable elimination 
via the restriction operator. In fact, since a restricted variable can only occur in its scope, its value can be 
determined when solving the subproblem it encloses. For instance 

{x\){{x2)A{xi,X2) II {x2,)B{xi,X2)) 
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represents a problem with three variables xi,X2,X3 and two atomic problems A,B, sharing xi. One can 
solve the subproblem {x2)A{x\,X2) w.r.t. X2 right away, thus eliminating X2 and obtaining a solution 
A'{x\) parametric w.r.t. xi; similarly for {x2,)B{x\,X'i)- Then x\ can be eliminated in (.ri)(A'(.ri) || 
B'{x\)), which yields the global solution. 

The specification provides axioms that are common in process algebras, and have a natural inter¬ 
pretation in terms of optimization problems. For instance: subproblems can be solved in any order, the 
formal names of their variables are irrelevant and so on. In particular, scope extension (actually, scope re¬ 
duction) becomes a crucial operation: it allows variables to be eliminated earlier, when solving a smaller 
subproblem. This produces a more efficient solution of the secondary optimization problem of dynamic 
programming. 

The computation of an optimal solution is formalized as an evaluation of terms in a suitable domain 
of cost functions, giving a cost to each assignment of free variables. This domain is indeed an algebra for 
the specification: constructors are interpreted as operations over cost functions, performing optimization 
steps, and axioms become useful properties. In particular, they tell that cost functions are preserved under 
rearranging the structure of the problem to get a more efficient computation. Moreover, the underlying 
nominal structure provides a notion of support, that is the set of variables that are really relevant for the 
computation. This is the key for a finite, efficient representation of cost functions. 

Then we introduce a graphical notation for optimization problems, called nominal hypergraphs. They 
are hypergraphs with an interface where names are assigned to some vertices. These vertices are vari¬ 
ables for the overall problem. Other vertices are variables of subproblems that must be eliminated in the 
optimization process. 

We show that nominal hypergraphs can be given an algebraic structure, and that (isomorphic) nomi¬ 
nal hypergraphs can be described by (congruent) terms of our specification. This allows us to recursively 
compute the cost function of a nominal hypergraph by performing the computation on any of the corre¬ 
sponding terms. Moreover, we show that well-known structures to represent parsing of graphs, namely 
tree decompositions Ealia, can be represented in our framework as terms. 

However, nominal graphs still lack a description of the variable elimination strategy. We describe 
this information as hierarchical nominal hypergraphs, that are trees describing the decomposition of a 
nominal hypergraph in terms of nested components, each corresponding to a subproblem. Such trees 
correspond to terms without the scope extension axioms, i.e., where the scope of restrictions is fixed. A 
bottom-up visit of a hierarchical nominal graph yields a dynamic programming algorithm with the given 
variable elimination strategy. 

We apply our approach to the e-mobility case study from the Autonomic Service-Component Ensem¬ 
bles (ASCENS) project. In ASCENS, systems are modeled as self-aware, self-adaptive and autonomic 
components running in ensembles (dynamic aggregations of components), which through interactions 
among them and with the environment accomplish both individual (local) and collective (global) goals 
by optimizing the use of resources. 

In the e-mobility case study, carried on in collaboration with Volkswagen, the traffic system is mod¬ 
eled as ensembles of electrical vehicles with the goal to optimize the usage of resources (electricity, 
parking places, etc.), while ensuring the fulhllment of individual goals (such as reaching in time the 
destination) and collective goals (avoiding traffic jams or guaranteeing that all vehicles hnd a spot where 
to park). So in general, besides optimizing local resources, for example by Ending the best trips and 
journeys for each vehicle, the e-mobility case study aims at solving global problems, involving large en¬ 
sembles of different vehicles. In ifTSl several optimization problems are presented for the e-mobility case 
study. In fS ] a parking problem is considered. In the formulation we present in this paper, parking sys¬ 
tems are regarded as nominal hypergraphs, where each hyperedge is a parking zone and vertices are cars 
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that may be parked inside the zone. Their term representation is interpreted as functions telling the cost 
of parking cars inside or outside certain parking zones. Thus cost functions only have binary arguments. 
While more efficient, this domain choice yields operators remarkably neater than those of the classical 
point-wise interpretation and our framework is flexible enough to accommodate them. 

2 An algebraic specification for optimization problems 

We introduce an algebraic specification for describing the structure of optimization problems. Variables 
are represented as names, belonging to an enumerable set M. We write Perm{M) for the set of permuta¬ 
tions over M, i.e., bijective functions n-M ^ M. 

Definition 1 (Optimization signature). Let C be a set of constants denoting atomic problems, equipped 
with an arity function ar: C ->• N telling how many variables each problem involves. We assume the empty 
problem nil, with ar{nil) = 0. The optimization signature is given by the following grammar 

p,q-= p II ^ I {x)p I pit I A{x) I nil 

where A^C, 7i € Perm{M), {x} ux<~J\f and |x| = ar{A) (we overload the notation x to indicate both a 
vector and a set of names). 

Here: 

• the parallel composition p \\ q represents the problem consisting of two subproblems p and q, 
possibly sharing some variables; 

• the restriction {x)p is p where the assignment for x has already been determined; 

• the permutation pit is p where variable names have been exchanged according to tt; 

• the atomic problem A{x) represents a problem that only involves the problem A over variables x; 

• nil represents the empty problem. 

We assume restriction has precedence over parallel composition. 

Free names of p are recursively defined as follows 

fn{p\\q)= fn{p)vjfn{q) fn{{x)p)=fn{p)\{x} fn{pn) = 7l{fn{p)) 

fn{A{x))=x fn{nil) = 0 

We consider syntax up to structural congruence axioms shown in fig.[T] The operator || forms a com¬ 
mutative monoid, meaning that problems in parallel can be solved in any order (AX||). Restrictions can 
be a-converted (AX^), i.e. names of assigned variables are irrelevant. Restrictions can also be swapped, 
i.e., assignments can happen in any order, and removed, whenever their scope is nil (AX(;f)). The scope 
of restricted variables can be narrowed to terms where they occur free (AXse)- Axioms regarding per¬ 
mutations say that identity and composition behave as expected (AX;;;) and that permutations distribute 
over syntactic operators (AX^). Permutations are assumed to behave in a capture avoiding way when 
applied to (x) p. We call optimization algebraic specification the specification made of the optimization 
signature and the congruence axioms, and optimization terms the terms for the specification. 

We include permutations in the specification because they provide a general mechanism to compute 
the set of “free” names in any algebra, called (minimal) support. 

Definition 2 (Support). Let A be an algebra for the optimization specification, and let be the inter¬ 
pretation of 71 in A. We say that X cj\f supports a e A whenever, for all permutations 7l acting as the 
identity on X, we have aTl^ = a. The minimal support supp{a) is the intersection of all sets supporting a. 

For instance, let be the interpretation of 7i on optimization terms: given a term p, ptf applies 7t to 
all free names of p in a capture avoiding way. It is easy to verify that supp{p) = fn{p). 
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(AX||) p\\qBq\\p {p\\q)\\r = p\\{q\\r) p\\nil = p 
(AX(^)) {x){y)p = {y){x)p {x)nil = nil 

(AXa) {x)p = {y)p[x<^y] {y i fn{p)) 

(AXse) {x){p II q) = {x)p II <7 {xi fn{q)) 

(AX;r) p\6 = p {pTi')n = p{TioTi') 

p A{xi,...,x„)K=A{K{xi),...,K{xn)) nilTl = nil {p\\q)K = pK\\qK 
^ {{x)p)n = {x)pn' {n'(x) = x,7t'{y) = n{y) for xi^y) 


Figure 1: Structural congruence axioms of the optimization specification. 

2.1 Hierarchical optimization specification 

The scope of restrictions determines a solution for the secondary optimization problem, because it spec¬ 
ifies when restricted variables should be eliminated. However, the presence of (AX^^) identifies ferms 
corresponding fo differenf solufions. We call hierarchical optimization specification fhe opfimizafion 
specification wifhouf (AXse), and hierarchical terms ifs freely generafed terms. 

We are interested in two forms of hierarchical terms. 

Definition 3 (Normal and canonical forms). A term is said to be in normal form whenever it is of the 
form 

(x)(Ai(xi) ||A2 (x2) II ••• \\An{Xn)) 

withAi&C(i= l,...,n) and x £ f j u • ■ ■ u It is in canonical form whenever it is obtained by the repeated 
application to a non-hierarchical term of {AXse), from left to right, until termination. For both forms, we 
assume that subterms of the form {x)nil (where x may be empty) are removed using (AX(_^^) and (AX||). 

Normal and canonical forms are somewhat dual: normal forms have all restrictions at the top level, 
whereas in canonical forms every restriction (x) is as close as possible to the atomic terms where x occurs 
(if any). A term in normal form is intuitively closer to a typical optimization problem: x specifies which 
variables should be assigned, and the term in its scope represents subproblems and their connections. 
In a term in canonical form, variables are eliminated as soon as possible. Notice that a term may have 
more than one canonical form, whereas normal forms are unique (up to the hierarchical optimization 
specification congruence). 

Remark 1. Hierarchical terms in normal and canonical form can be regarded as canonical represen¬ 
tatives of =-classes (recall that = is the structural congruence of fig. [^, because = is coarser than the 
hierarchical optimization specification congruence. 

3 Optimization problems as nominal hypergraphs 

Recall that a hypergraph G is a triple {Vg,£ 0 , 00 -Eg V^), where Vg is the set of vertices. Eg is the set 
of hyperedges and, for each e e Eg, aG{e) is the tuple of vertices attached to e (F^ is the set of tuples over 
Vg)- Let She. a set of edge labels, equipped with a function ->• N telling the number of vertices ar{l) 
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of an edge with label 1 . A labeled hypergraph G is a hypergraph G plus a function lab-E g £ mapping 
each hyperedge e € Eg to its label I such that \aG{e) \ = arl{l). Given two (labeled) hypergraphs Gi and 
G2, we write Gi u G2 for their component-wise disjoint union. 

Optimization problems can naturally be seen as hypergraphs labeled over atomic subproblems, where 
vertices correspond to variables. We introduce a notion of labeled hypergraph where some vertices are 
associated variable names. 

Definition 4 (Labeled nominal hypergraph and their morphisms). A labeled nominal hypergraph (NEl- 
graph in short) is a pair TfJ > G, where G is a labeled hypergraph without isolated vertices and rj is a 
partial injection from Vg to M, assigning names to some vertices of G. The set img{r\) is called the 
interface of G and def{ri) (the domain of definition of rj) are called interface vertices. Given two NH- 
graphs T]i > G\ and r\2 t> G2, a NH-graph morphism h- T]i > Gi ^ T]2 > G2 is a homomorphism Gi ->• G2 of 
labeled hypergraphs that preserves names, namely r\\ohy = T]2, where hy is the action ofh on vertices. 

Interface vertices can be understood as “external” vertices, with a public, global identity. They may 
be interaction points, i.e., they may be shared, with other graphs. This will allow for a simple definition of 
parallel composition of NH-graphs. Notice that NH-graph homomorphisms must be injective on vertices, 
because they must commute with functions that are injective on vertices. 

We say that rp > Gi and T]2 > G2 are isomorphic, written rp > Gi = r]2 > G2, whenever there is an 
NH-graph isomorphism (i.e., a NH-graph morphism whose underlying hypergraph homomorphism is an 
isomorphism) between them. 

Remark 2. A NH-graph rj > G can be seen as the following span of (total) injective graph homomor¬ 
phisms 

m^^[img(r])f^G 

where [img(ri)] is the discrete graph with vertices img{r\), J\f is the infinite discrete graph with vertices 
J\f, t]/(v) = T](v) and rjr A an embedding. 

3.1 Example 

Consider the optimization term shown in the introduction, without the outer restriction 

(X2)A(X1,X2) II {x 2 ,)B{xi,X 2 ) 

Recall that such a term may represent the following optimization problem: given two subproblems A and 
B, with cost functions parametric in x\, X2, and x\, X3, respectively, find the optimal total cost. Actually, 
here xi is free, meaning that the total cost is parametric in xi. 

The problem can be represented as a NH-graph rj >G with two hyperedges, labeled A and B, and three 
vertices vi,V2,V3, corresponding to xi,X2 and X3. Actually, only xi becomes an 
interface name, because it is the only variable the problem “exposes”. Other vari¬ 
ables are not part of the interface, meaning that they are taken up to a-conversion. 

The NH-graph is depicted on the right: the dashed line describes the domain of 
definition of rj , namely T] (vi) = xi. 

3.2 An algebra for NH-graphs 

Now we show that we can regard NH-graphs as elements of an algebra for the optimization specification. 
This will allow us to recursively evaluate (parsing of) NH-graphs as cost functions, as done in section]^ 
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Operations are interpreted as follows 


T]l > Gi P T72 > G2 = T] 0 (Gi '±!G 2 )l^y 


r^y = 


|(vi,V2) e Vgi Vgj 


T?l(vi) = 772 (v2) 
t undefined 


where • 




|t]i(v) V€Vg, 
[ t ] 2 ( v ) V€Vg2 


{xy{ri > G) = Tjxx 0 G where 
( t ] > G)ti^ = {7iori)>G 


( undefined ri{v)=x 
ri(v) ofherwise 


The parallel composifion T]i > Gi ||g T]2 > G2 is compufed by faking fhe disjoinf union of fhe fwo NH- 
graphs and fhen idenfifying vertices wifh fhe same inferface names (formally, / fakes equivalence 
classes of vertices). The funclion 77 is defined on equivalence classes of verfices as expecfed. The restric¬ 
tion (x)gT] > G of T] > G w.r.t. X simply removes x from the interface of T] > G. 

The interpretation of constants can be defined via a mapping 


A 

lA{xi,X2,...,Xn)f = ... Inilf = l- 0 -^J\ft>OG 

V V V 

Xi X2 Xji 

where !: 0 ^ AA is fhe (unique) mapping from 0 to J\f and Og is fhe empfy hypergraph 

Now we have fo show fhaf we have indeed defined an algebra. We have fo check fhaf congruence 
axioms are safisfied. We firsl need fhe following characferizafion of fhe minimal supporf of a NH-graph. 

Lemma 1 . supp{r] t>G) = img(ri). 

Proposition 1 . Operations P, (x)^, 7 t^ satisfy axioms of fig. where = becomes = and fn{-) becomes 

supp{-). 

Now we can define a unique evaluation of opfimizafion terms: given p, fhe corresponding NH-graph 
can be compufed by sfrucfural recursion using fhe evaluafion of constanfs given above and fhe infer- 
prefafion of operafions on NH-graphs. This induces a sound and complefe axiomafizafion for NH-graphs. 
In facl, sfrucfurally equivalenf opfimizafion ferms are evaluafed fo isomorphic NH-graphs (soundness). 

For completeness, given a NH-graph rj > G, we can construcl an equivalenf ferm in normal form 
(regarded as canonical represenfafive of ifs E-class, by remark[T]) as follows. We encode each edge e of G 
as an atomic ferm wifh label labG{e). The argumenfs of fhis ferm are fhe names of fhe inferface verfices 
of e, and arbifrary ones for non-interface verfices. Finally, we form fhe parallel composifion of all fhese 
terms and we resfricf fhe names fhaf are nol in fhe inferface. Every choice of resfricfed names is valid: 
a-conversion guaranfees fhaf all possible encodings of T] > G are in fhe same sfrucfural congruence class. 

3.3 Tree decompositions 

In graph theory, we have the well-known notion of tree decomposition of a graph ll 22 l . which can be 
understood as a way of parsing a graph. We report the definition by ifT^ . 
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Definition 5 (Tree decomposition). A tree decomposition (TD) of a hypergraph G is a pair {T,X), where 
T - {N,A) is a tree (i.e., an undirected, acyclic graph) and X = {Xn}mN i^ a family of subsets X„ c Vg, 
one for each node of N, suchthat: (a) U«eA?^n = VG; (b) for all e ^ Eq there is n ^ N such that aciE) ^X,,; 
(c)for all V € Vg, the set of nodes {n€N\v€ X^} induces a subtree ofT. 

In our framework, ordinary hypergraphs (without isolated vertices) can naturally be seen as NH- 
graphs with empty interface, and each of their TDs can be represented as an optimization term. 
Theorem 1. Every TD for G induces an optimization term t such that = t > G, where V-Vg^ ff is 
the nowhere defined function. 

The idea is constructing t via a visit of T from a chosen root r. We first associate the induced subgraph 
G„ of G to X„, for each node n € N. Each time a new node n is expanded in the visit, we generate 
the following subterms of t, all in parallel: one term representing edges and nodes of not already 
in f, the subterms corresponding to n’s children. Correctness of the translation is guaranteed by (a)- 
(c) of definition Notice that, since all variables are restricted, choosing a different root amounts to 
rearranging restrictions. By soundness, this operation results in terms with isomorphic images via 
all isomorphic to G. 

4 Representing and solving optimization problems 

We now show how typical optimization problems can be represented and solved in our algebraic frame- 
workj^Suppose we have n atomic problems Ai ,... ,A„ whose variables can be assigned values in D, and 
we want to minimize a function of the form 

E 

\<i<n 

where each ^ Moo, for / = 1,... ,n, gives a cost to each variable assignment for the problem 

A, ; an infinite cost represents a forbidden assignment. 

The problem can be represented as the following term in normal form 

p = {x){Ai{xi) \\ ■■■ \\ A„{x„)) where x = xiU---ux„ 

and the computation of the optimal cost as a function 

[[p]]'^:(Ar ->• D) Moo 

giving a cost to each assignment of variables. More precisely, its computation is performed by assigning 
values to the free variables of p (discarding assignments to other variables), and minimizing w.r.t. bound 
ones. Typically we have fn{p) = 0, so minimization is performed w.r.t. all variables. 

Formally, we take an algebra for the optimization specification formed by cost functions (j ): {M 
D) ->• Moo, where we interpret optimization terms. For any assignment of variables p:A/^ D, the inter¬ 
pretation of constants is 

lAi{xi,...,Xn)f p =CA,{p{xi),...,p{x„)) lnilfp = 0 

and complex terms are recursively interpreted as follows 

bi II Pifp = IPifp + lPifp Kx)pfp^mmlprip[x^v]) IpTtfp = IpfipoK-^) 

veB 

We have the following property, which comes from the theory of permutation algebras. 

* A different, more efficient, setting will be described in section 5.1 
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Property 1. supp{lpf)cfn{p). 

We introduce a condition on cost functions, called compactness. A compact (j) depends only on a 
“few” variables. This is essential to compute and store cost functions in an efficient way, as we will see 
later. 

Property 2 (Compactness). We say that (j)-.{Af ->■ D) ^ Mc>o A compact if P\supp{(j>) = P'\supp(<i)) ifnplies 
0 p = ^p',foranp,p'.M^B. 

When considering a term p, by property]^ is compact if it depends only on assignments to free 
variables of p. This is clearly the case for the interpretation of constants, and can be shown by structural 
induction for complex terms. Notice that this property is not true for the whole algebra of functions 
{M ^ D) ^ ]Roo^ but only for the subalgebra in the image of 
Canonical forms and normal forms of a term always have the same cost function. This is thanks 
to the following proposition, which is a direct consequence of cost functions forming an algebra of the 
optimization specification. 

Proposition 2. If p = q then \pf = [[^]]‘^. 

4.1 Computational complexity of cost functions 

Although structurally congruent terms have the same cost functions, these functions may be computed 
in different ways, each possibly with a different computational cost. In fact, the position of restrictions 
inside a term determines a strategy for variable elimination. As already mentioned, finding fhe besf one 
amounfs fo giving a solufion for fhe secondary opfimizafion problem. 

We infroduce a nofion of complexify for an opfimizafion problem p, similar fo fhe one of |Tl, esfimaf- 
ing fhe cosf of compufing ifs value This is given by fhe funcfion wifh greafesf “size” encountered 
while inducfively consfrucfing [pj^, fhe size being given by fhe variables in fhe support, fhaf are fhe only 
ones defermining fhe value of IpY (propertyj^. 

Formally, fhe complexify of p, wriffen ((p)), is recursively defined as follows 

((A(r))) = |r| {{nil)) = 0 (((r)p)} = ((p)) ((p || ^)) = max{((p)), ((^)),|/n(p || ^)|} 

The inferesfing cases are {x)p and p || q\ fhe compufafion of |I(r)p]]‘^ relies on fhaf of [pj^, whose support 
may be bigger, so we sef fhe complexify of (r)p fo fhaf of p; compufing \p || q}f requires compufing [pj^ 
and buf fhe supporf of fhe resulfing funcfion is fhe union of fhose of p and q, so we have fo find fhe 
maximum value among ((p)), ((^)) and fhe overall number of free variables. 

Complexify is well-defined only for hierarchical terms: applying (AXse) to choose a different vari¬ 
able elimination strategy may change the complexity. Consider, for instance, the following term in normal 
form 

p={x\){x2){x-i){A{xi,X2) \\B{x2,x-i)) ; 

we have ((p)) = 3 , but if we take a canonical form 

q={x2){{xi)A{x]_,X2) II {xf)B{x2,xf)) 

we have ((^)) = 2 . Indeed, we have the following results for hierarchical terms. 

Lemma 2. Given (x)(p || q), with x ^ fn{q), we have {{{x)p || q)) < (((r)(p || q))). 

As an immediate consequence, all the canonical forms of a term always have lower or equal com¬ 
plexity than the normal form. 

Theorem 2. Given a term p, let n be its normal form. Then, for all canonical forms c of p we have 
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(a) |A(;ci,;c 2 )r 


(b) lB{x2,X3)f (c) |( 2 ci)A(;ci,X 2 )f 


X\ X2 cost 

d\ di 1 

d\ d2 5 

d 2 d\ oo 

d2 d2 2 


X2 X2 cost 

d\ d\ 9 

d\ d2 1 

d2 d\ 6 

d2 d2 13 


X2 cost 

d\ min{ 7 ,oo }=7 
d2 min{ 5 , 2 } = 2 


(d) |(x3)B(x2,X3)f 


(e) |(xi)A(xi,X 2 ) II {x 2 ,)B{x 2 ,x-i)f 


X2 cost X2 cost 

di min{ 9 ,l} = l d\ 7 + 1 = 8 

d2 min{ 6 , 13 } = 6 d2 2 + 6 = 8 


Table 1 : Cost functions for the problems in the example. 


4.2 Example 

Consider two problems A and B, with two variables each, ranging over {d\,d2}- Their cost functions are 
shown in Tables [T^ and We consider the optimization problem that consists in finding the minimal 
value of A(xi,X2) +B(x2,X3). 

As we already saw, the term in canonical form representing the problem is 


p = {x2){{x\)A{xi,X2) II {x3)B{x2,X3)) 


We now show how can be computed. We proceed in a bottom-up order, from atomic subterms to 
increasingly complex terms. This is close to a dynamic programming algorithm, as it allows computing 
and storing a (finite, thanks to the compactness property) representation of cost functions once and for 
all. Table [T] show such finite representations in a tabular form. We perform the following optimization 
steps, each corresponding to an operator of the syntax: 

1. \{xi)A{x\,X2)f and \{x 2 )B{x 2 ,X 2 ,)Y computed by minimizing [[A(xi,X2)]]'^ and \B{x 2 ,X 2 )Y 
w.r.t. x\ and X3 respectively (Tables and [Td| ). Notice that these functions can be computed in 
parallel. 

2. \{xy)A{x\,X2) II {x2)B{x2-,xj,)Y is computed by evaluating \{x\)A{x\,X 2 )Y and \{xj,)B{x2-,X'i)Y 
on the same value for X2, and then summing up the results (Table [Te]). 

3 . Finally, \{x2){{x\)A{x\,X2) || {x3)B{x2,xj))Y is computed by minimizing the function of step 2 ) 
w.r.t. X2- 


The last step gives the overall minimal value 8. By looking at Tables in a top-down order, from le 


to each time picking those variable assignments that contributed to the cost, one can recover the 
corresponding optimal assignment(s) for xi,X2 andx3, namely d\,d\,d2 and d2,d2,d\. 


5 Dynamic programming on hierarchical NH-graphs 

The existence of an algebra of NH-graphs allows us to recursively compute cost functions for these 
graphs: given a NH-graph rj > G and an optimization term such that - rj > G, we can compute its 
cost function in the style of section 
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{x2){{xi)A{xi,X2) II {x3)B{x2,X3)) 


0 


© 





Figure 2 : A hierarchical term and the corresponding hierarchical NH-graph. 


However, the information about the variable elimination strategy cannot be recovered from the NH- 
graph itself. In order to do this, we need to introduce the graphical counterpart of hierarchical terms, 
which we call hierarchical NH-graphs. They are trees that describe the structure of a NH-graph rj >Gm 
terms of nested components. These trees are such that: 

• the root is the discrete hypergraph formed by the interface vertices of G; 

• each internal node n is a discrete subgraph of G; 

• leaves are hypergraphs with a single hyperedge of G; 

• there is an arc from G to G' whenever G c G'. 

The intuition is that each internal node n of the tree is a component of t] > G that exposes some additional 
vertices and includes all the components in the subtree rooted in n. Leaves are basic components, i.e., 
hyperedges. 

The correspondence between hierarchical terms and hierarchical NH-graph graphs is exemplified in 
fig-@ The scope of each restriction determines a component in the tree, where a vertex for the restricted 
name is added. For convenience, we used the same name for restricted variables and corresponding non¬ 
interface vertices, but the latter, as in ordinary NH-graphs, are actually up to a-conversion. A top-down 
visit of the tree amounts to “opening” scopes and revealing their names. 

As hierarchical terms, hierarchical NH-graphs describe a solution for the secondary optimization 
problem. It is possible to show that hierarchical NH-graphs form an algebra for the hierarchical op¬ 
timization specification. As seen in the example, each hierarchical term corresponds to a hierarchical 
NH-graph. We also have the opposite correspondence. 

Proposition 3 . Each hierarchical NH-graph can be represented as a hierarchical term. 

Exploiting this correspondence, we can interpret hierarchical NH-graphs as evaluations of cost func¬ 
tions with a specific variable eliminafion sfrafegy, which can be implemenfed via dynamic programming. 

Remark 3 . In remark we have characterized NH-graphs as spans of hypergraph homomorphisms. 
Interestingly, exploiting this characterization we can recover the NH-graph of which a given hierarchical 
NH-graph is the decomposition. In fact, a hierarchical NH-graph can be regarded as a span where the 
right part is a diagram T made of a “tree” of graph embeddings. 


m^[img{r])Y^T 
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Here rjr maps interface vertices to themselves in the root of T. Notice that the morphisms of T tell which 
are the interface vertices in the nodes ofT. In order to “paste” together leaf hypergraphs ofT, we can 
make a colimit of T in the category of graphs and their morphisms. The result is the disjoint union of 
such hypergraphs, where vertices that are images of the same interface one, along the morphisms of T, 
are identified. 

5.1 The parking optimization problem 

We now introduce the parking optimization problem and we apply our approach to it. The parking op¬ 
timization problem consists in finding the best parking zone for each vehicle of an ensemble, that is a 
group of vehicles with similar features. This can be formalized as follows. Assume a set of parking zones 
P = {A,B,...} and of car variables C = {x,y,...}, and two functions: 

• c:P N, assigning a capacity to every zone; 

• F'.C ^ P ^ Moo, specifying the cost F{x){A) for x to park in A. 

Given an assignment p:C ->• P of cars to zones, let Pa = {x\ p{x) = A}. We want to find an assignmenf p 
such fhaf IpaI < c(A), for all A e P, minimizing 

E F{x){p{x)) 

xeAf 

The problem can be specified in fhe sfyle of section Here a ferm p represenfs a parking sysfem: 
A{x\ ,...,x„) means fhaf jc, mighf be parked in A; (x) p means thaf car a: cannof be parked oufside of p, so 
if musf have a parking spof in one of fhe zones of p. In general, a ferm p represenfs a parf of fhe sysfem 
made of one or more parking zones. 

In section]^ we presenfed an algebra of cosf functions for fypical opfimizafion problems. Here we 
infroduce anofher, “more efficienl” algebra which, neverfheless, fils info our algebraic framework and 
can be used to evaluate optimization terms. To each parking system p we associate a function 

IpT ■'Pifn{p)) 

The intended meaning of is the cost of parking in p cars X c fn{p) and all cars corresponding to 

variables restricted in p. Notice that the evaluation function is quite different from that of section|^ 
Here assignments do not fix fhe values of variables, i.e., fhe parking zones where cars are allocafed, buf 
only fheir positions wifh respecf fo fhe presenf p. However, our framework is able fo accommodafe also 
fhis more efficienl selling. 

To avoid handling polymorphic funclions, we aulomafically exfend funcfions \pY to the whole M, 
namely \pf'-P{-^) by letting 

lpfX = lpf{Xnfn{p)). ( 1 ) 

However, this function is still determined by subsets of fn{p), so it admits a finite representation that 
can be efficiently computed and stored. Formally, properties [T] and hold for extensions to M (regarded 
as functions {J\f ->-{ 0 , 1 })^ Moo). 

Cost functions can be defined by recursion on fhe slrucfure of syslems. As menlioned, if is enough fo 
define fheir acfion on a subsel X of fheir supporl. We have 


|X|<C(A) 

lA(x)l X = \ , 

oo otherwise 
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meaning that cars X can be parked in zone A iff their number is at most the capacity of A. For nil we 
simply have = 0. 

Then we have 


Ip II qfX = min 

{XuX2hp2(X) 


lpfX, + lqrX2 


Xi ^fn{p), 1 

X 2 ^fn{q) J 


where V2{X) are the partitions in two sets of X c fn{p || q). Here, to park cars X in component p || q, 
one has to park each of them in either component p or component q, but not in both. Thus the best option 
must be chosen. Finally, we have 

\{x)pfX=lpf{^^{x})- 

Here it is required that car is parked in component p. 

Typically, the whole system s has no free names. Thus |[5]|0 is a real number, the total minimized 
cost, or oo if the problem has no solution. 

In order to have a proper theory of cost functions, we have to show that we have indeed defined 
a model of the optimization specification. Let Xn he. the element-wise application of a permutation 
71: A/^ ->• AT to X c J\f. Then we have the following theorem. 

Theorem 3. Cost functions <f'-V{M) Mcxi satisfying (Q/orm a model of the optimization specification, 
together with the given interpretation of operators and the permutation action {(j) 7 t)X = (j){X 7 l^^). 


5.2 Dynamic programmin g algorithm 

Consider the scenario with three possible parking zones A,B,C and three cars xi,X2 and V3. We assume 
the following values for F and c. 


F(vi)(A) = 3 
^’(x2 )(A) = 4 
F{x3){A) = 00 


F{xi){B) = 00 
F{x2){B) = 6 
F{x3){B) = 4 


F{xi){C) = 00 
F{x2){C) = 00 
F{x3){C) = 1 


c(A )=2 c{B) = 2 c(C) = 2 


In fig. on fhe leff side, we show fhe ferm in normal form, and fhe corresponding NH-graph, mod¬ 
eling fhe system. We wanf fo compufe fhe cosf function \pY using dynamic programming. The crucial 
properly is compacfness of cosf functions: \pY can be represented as a fable of size |/n(p)|^ Allhough 
a problem is fypically specified as a term in normal form, we consider ils canonical form c, shown on fhe 
righf side of fig.[^ because ils complexify is equal or lower (Iheoremj^. We assume lhal all occurrences 
of nil have been eliminated via slruclural congruence. 

We propose a dynamic programming algorilhm lhal is driven by fhe hierarchical NH-graph lhal c 
describes, shown in Ihe righl side of fig. This algorilhm operates as fhe one in section [ 4 ~ 2 l bul fables 


are computed using fhe differenl inlerprelalion of operators on cosl functions, inlroduced in section 5.1 


The algorilhm slarls from Ihe cosl functions for Ihe leaves. These are shown in Tables [ 2 a|to|^ where 
Ihe leflmosl columns indicates whelher a car is parked inside (-/) or oulside (-) each zone. They are 
computed as described in section 5.1 e.g.,for |[A(xi,X2)]|^, the cosl for each row is |[A(xi,X2)]|‘^2f, where 
X are Ihe variables marked wilh •/ in lhal row. Then, Ihe algorilhm performs a bollom-up visil of Ihe free 
and eliminates variables accordingly. More precisely, whenever an edge from G to G' is Iraversed, wilh 
G and G' discrete hypergraphs, variables G \ G' are eliminated. In Ihe following we show Ihe elimination 
steps, also indicated in fig.|^ 

^The actual size is gjiow the exponent, as 2'* < 2^ iff x< y. 
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Normal 


p = {xi){x2)ix3){A{xi,X2) \\B{x2,X3) || C(;C3)) 


Canonical 

C={x2){{xi)A{x\,X2) II {x3){B{x2,X3) || C(;C3))) 


0 


A 

B 


C 

Ms 

/ 




hi 


42) 


© © 


© © 



B 

© © 


c 


Figure 3 : Graphical and corresponding term representation of a parking problem. 


( 1 ) Elimination of xp. Table \{x\)A{x\,X2)\^ 
to be inside A; 


, with only one column X2, is computed by forcing x\ 


( 2 ) Elimination of X3\ the table \{x3){B{x2,X3) || C(.r3))]|^ is computed: Table 2 e shows values for ;c2. 
the partitions considered when computing the output cost, and the final cost. Notice that this and the 
previous step could be executed in parallel. This fact comes immediately from terms {x\)A{x\,X2) 
and {x3){B{x2,X3) || C{x3)) being composed in parallel in {x\)A{x\,X2) || {x3){B{x2,X3) || C{x3)). 

( 3 ) Elimination of xg- finally, the Table \pY ( pf| ) is computed, by comparing costs of parking X2 inside 
either (xi)A(xi,X2) or {x3){B{x2,X3) || C{x3)). 


( 4 ) Optimal variable assignment: tracking back through the Tables we find: 

• X2 inside {xi)A{xi,X2) || {x3){B{x2,X3) || C(x3)); 

• X2 inside {x\)A{x\,X2)', 

• X2 inside A{x\ ,^2) with cost 4 ; 

• xi inside A(xi,X2) with cost 3 ; 

• X3 inside B(x2,X3) || C{x3)\ 

• X3 inside C{x3) with cost 1. 


Notice that, in general, the outcome of the algorithm may be 00, whenever there is no car assignment to 
parking zones that respect capacities. 


6 Conclusion 

In the paper we have introduced two process algebra-like specifications for the description of optimiza¬ 
tion problems. The more abstract version (which includes the scope extension axiom) defines (hyper) 
graphs, where vertices are variables, and edges are tasks whose costs depend on the values of the adja¬ 
cent variables. Dropping the above axiom yields a specification corresponding to different parsing trees 
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Table 2 : Example tables. Parameters of atomic subterms are often omitted. 


(a)| 

lA(xi 

A2)r 

(b) 

{B{x 2 

A3)r 

(c) [c(x 3 )r 

Xl 

X2 

cost 

X2 


cost 

2C3 

cost 

v/ 

xA 

1 

sT 

xA 

10 

x/ 

1 


- 

3 

xA 

- 

4 

- 

0 

- 


4 

- 

xA 

6 



— 

— 

0 

— 

— 

0 




(d) |(xi)A(xi,X2)r (e) |(x3)(B(x 2,X3) || C(x3))f 


X2 

cost 

2C2 

m 

C 

icf 


cost 

s/ 

1 


X3 

Xl 

2C3 



— 

3 



xA 

- 

10 




xA 

- 

xA 

xA 

7 

1 




xA 

- 

- 

4 

1 




- 

- 

xA 

1 

1 


(f) I('*2)( (xi)A(xi,X2) II (x3)(B(x 2,X3) || C(x3) ) )f 


I(M)Ar 

II c)r 

[(xi)Af- i-|[(x 3)(B II C)f cost 

X2 

X2 


xA 

- 

8 

- 

xA 

30 

o 


of the given graph. Choosing a particular tree corresponds to selecting a dynamic programming strategy 
for the given problem, whose execution can be carried on via a bottom up visit of the tree. We apply 
our approach to the parking optimization problem developed, in collaboration with Volkswagen, in the 
ASCENS e-mobility case study. 

The idea of exploiting graphs to decompose and solve various kinds of problems is not new. In ifTTI 
graphs are represented as elements of an algebra and monadic second-order properties are evaluated on 
them. In Q dynamic programming algorithms are derived from (nice) tree decompositions of graphs. In 
['51 tree decompositions are represented as in a category of spans and cospans, and then as terms of an 
algebraic specification. Our approach has the following advantages w.r.t. the cited ones: 

• Our algebraic specification is simpler, but nonetheless expressive. In fact, variable elimination 
strategies can be represented via restrictions. Moreover, we have a graphical representation of such 
strategies as hierarchical NH-graphs, which can be regarded as very simple tree decompositions. 

• Our algebras are permutation algebras, which provide: (a) a state-of-the-art treatment of a- 
conversion and of freshness requirements; (b) a uniform and general definition of domain given by 
the notion of support. Operations are defined on the whole set of names, so they are independent 
of the actual interface (support), unlike ||6l. Moreover, the notion of support automatically defines 
the sizes of the tables employed in the dynamic programming implementation. 

The following lines of research are also related to our work. Bistarelli, Montanari and Rossi deal 
with SCSP | 2 i and its combination with logic programming IS) lU and concurrency El ■ They give an 
interpretation of constraints over certain semirings, e.g., the tropical semiring, as we do here. However, 

























94 


Dynamic Programming on Nominal Graphs 


operations on constraints are defined point-wise using the semiring operations: the approach is too re¬ 
strictive, e.g., it does not easily accommodate the case study shown in this paper. A direct connection 
between (logical) CSP and dynamic programming is shown in || 20 )I . Dechter in |[T 3 l introduces bucket 
elimination as a general solution technique for a variety of problems: it consists in a strategy of prob¬ 
lem reduction employing a convenient elimination ordering of variables and constraints. The associated 
technique of conditioning search allows for approximated versions of the bucket elimination approach. 
Kohlas and Pouly in ifTTl suggest valuation algebras as a foundation for a general view of information 
processing. They define axioms for valuation algebras, consider a number of instantiations and provide 
generic inference algorithms for their processing. Our approach is similar, but more direct, being based 
on a simple process algebra specification and on a bottom up visit of a tree of graphs satisfying the 
specification. In lITOl distributed systems are represented as CHARMs (Concurrency and Hiding in an 
Abstract Rewriting Machine), that are hypergraphs with a global and a local part. They form an algebra 
including edges and vertices restriction. Our algebras and NH-graphs are similar, but we do not need 
edges restriction. 

Dynamic programming evaluations are particularly interesting for autonomic systems, as studied by 
the ASCENS project, where the actual behavior often consists, typically for the dynamic programming 
case, of propagating local knowledge to obtain global knowledge and getting it back for local decisions. 
When dealing with global problems, however, the complexity of the dynamic programming algorithms 
can grow exponentially even for graphs of limited complexity. Consider a rectangular grid of size n, 
with vertices labeled by variables, and edges by cost functions with two arguments. It is shown in ifTSl 
that its complexity is exponential in n. There are efficient algorithms for finding the optimal elimination 
order of vertices in a graph, but they deal with specific cases (e.g., Gaussian elimination |[T 2 l| 23 l). Thus 
approximation techniques are quite relevant, in particular when a good global solution, possibly not 
optimal, is still acceptable. 

Several heuristic techniques can be experimented. For instance, for the parking problem we could 
restrict the number of possible zones for each car, taking the best k of them. Then if the optimal solution 
would include a choice worse than k for some car, the solution found, if any, would not be optimal. 
However, at least no client would be treated too badly. Another, quite general, approximation technique 
would be to artificially reduce the dimensions of tables by decomposing high dimensional ones into 
the sums of a few lower dimensional tables. The latter can be computed minimizing the mean square 
error ii 3 - The storage reduction can be propagated in such a way to reduce substantially the overall 
complexity. 

An interesting piece of future work would be to extend our approach to graphs which are incremen¬ 
tally modified, e.g., extended, at run time. The resulting scenario could consist of a (soft) (concurrent) 
constraint component together with a mobile pi-calculus-like process algebra component. A good exam¬ 
ple of this combination is cc-pi f 91 - Other aspects should be investigated in a precise way: the correspon¬ 
dence between classes of terms (normal,canonical) and NH-graphs; the nominal structure of NH-graphs. 


References 

[1] Umberto Bertele & Francesco Brioschi (1973): On Non-serial Dynamic Programming. J. Comb. Theory, Ser. 
A 14(2), pp. 137-148, doi: 10.1016/0097-3165(73)90016-2 

[2] Stefano Bistarelli, Ugo Montanari & Francesca Rossi (1997): Semiring-based constraint satisfaction and 
optimization. J. ACM44(2), pp. 201-236, doi: 10.1145/256303.256306 


N. Hoch & U. Montanan & M. Sammartino 


95 


[3] Stefano Bistarelli, Ugo Montanari & Francesca Rossi (2001): Semiring-based contstraint logic programming: 
syntax and semantics. ACM Trans. Program. Lang. Syst. 23(1), pp. 1-29, doi: 10.1145/383721.383725 

[4] Stefano Bistarelli, Ugo Montanari & Francesca Rossi (2002): Soft Constraint Logic Programming and Gen¬ 
eralized Shortest Path Problems. J. Heuristics 8(1), pp. 25^1, doi: 10.1023/A:1013609600697 

[5] Stefano Bistarelli, Ugo Montanari & Francesca Rossi (2006): Soft concurrent constraint programming. ACM 
Trans. Comput. Log. 7(3), pp. 563-589, doi: 10.1145/1149114.1149118 

[6] Christoph Blume, H. J. Sander Bruggink, Martin Friedrich & Barbara Konig (2013): Treewidth, pathwidth 
and cospan decompositions with applications to graph-accepting tree automata. J. Vis. Lang. Comput. 24(3), 
pp. 192-206, doi: 10.1016/j.jvlc.2012.10.002 

[7] Hans L. Bodlaender & Arie M. C. A. Koster (2008): Combinatorial Optimization on Graphs of Bounded 
Treewidth. Comput. J. 51(3), pp. 255-269, doi: 10.1093/comjnl/bxm037 

[8] Tomas Bures, Rocco De Nicola, Dias Gerostathopoulos, Nicklas Hoch, Michal Kit, Nora Koch, Gia- 
coma Valentina Monreale, Ugo Montanari, Rosario Pugliese, Nikola B. Serbedzija, Martin Wirsing & Franco 
Zambonelli (2013): A Life Cycle for the Development of Autonomic Systems: The E-mobility Showcase. In: 
SASOW, pp. 71-76, doi: 10.1109/SASOW.2013.23 

[9] Maria Grazia Buscemi & Ugo Montanari (2007): CC-Pi: A Constraint-Based Language for Specifying Ser¬ 
vice Level Agreements. In: ESOP, pp. 18-32, doi: 10.1007/978-3-540-71316-6_3 

[10] Andrea Corradini, Ugo Montanari & Francesca Rossi (1994): An Abstract Machine for Concurrent Modular 
Systems: CHARM. Theor. Comput. Sci. 122(1&2), pp. 165-200, doi: 10.1016/0304-3975(94)90206-2, 

[11] Bruno Courcelle & Mohamed Mosbah (1993): Monadic Second-Order Evaluations on Tree-Decomposable 
Graphs. Theor. Comput. Sci. 109(1&2), pp. 49-82, doi: 10.1016/0304-3975(93)90064-7 

[12] Elias Dahlhaus (2002): Minimal elimination ordering for graphs of bounded degree. Discrete Applied Math¬ 
ematics 116(1-2), pp. 127-143, doi:10.1016/S0166-218X(00)00331-0 

[13] Rina Dechter (1999): Bucket Elimination: A Unifying Framework for Reasoning. Artif. Intell. 113(1-2), pp. 
41-85, doi: 10.1016/S0004-3702(99)00059-4 

[14] Fabio Gadducci, Marino Miculan & Ugo Montanari (2006): About permutation algebras, (pre)sheaves and 
named sets. Higher-Order and Symbolic Computation 19(2-3), pp. 283-304, doi: 10.1007/sl0990-006-8749- 
3 

[15] Nicklas Hoch, Kevin Zemmer, Bernd Werther & Roland Siegwart (2012): Electric vehicle travel 
optimization-customer satisfaction despite resource constraints. In: IEEE IVS, pp. 172-177, 
doi: 10.1109/IVS.2012.6232240 

[16] Ton Kloks (1994): Treewidth, Computations and Approximations. Lecture Notes in Computer Science 842, 
Springer, doi: 10.1007/BFb0045375 

[17] Jiirg Kohlas & Marc Pouly (2011): Generic Inference: A Unifying Theory for Automated Reasoning. John 
Wiley & Sons, Inc., doi:10.1002/9781118010877.ch2, 

[18] Alberto Martelli & Ugo Montanari (1972): Nonserial Dynamic Programming: On the Optimal Strategy of 
Variable Elimination for the Rectangular Lattice. J. Math. Anal. Appl. 40, pp. 226-242, doi: 10.1016/0022- 
247X(72)90046-7 

[19] Ugo Montanari (1971): On the Optimal Approximation of Discrete Functions with Low-dimensional Tables. 
In: lEIP Congress (2), pp. 1363-1368. 

[20] Ugo Montanari & Francesca Rossi (1991): Constraint Relaxation may be Perfect. Artif. Intell. 48(2), pp. 
143-170, doi: 10.1016/0004-3702(91)90059-S 

[21] A. M. Pitts (2013): Nominal Sets: Names and Symmetry in Computer Science. Cambridge Tracts in Theoret¬ 
ical Computer Science 57, Cambridge University Press, doi: 10.1017/CB09781139084673 

[22] Neil Robertson & Paul D. Seymour (1984): Graph minors. III. Planar tree-width. J. Comb. Theory, Ser. B 
36(1), pp. 49-64, doi: 10.1016/0095-8956(84)90013-3 


96 


Dynamic Programming on Nominal Graphs 


[23] Mihalis Yannakakis (1981); Computing the Minimum Fill-In is NP-Complete. SIAM Journal on Algebraic 
Discrete Methods 2(1), pp. 77-79, doi; 10.1137/0602010 


